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The spectrum for the pure gauge sector is calculated for an extended Isgur-Paton model in 2+1 and 3+1 
dimensions and compared to recent lattice calculations of the glueball spectrum. The IP model is extended by 
inclusion of a rigidity (curvature) term and, in D=2+l, mixing through a higer topological contribution. For 
a choice of parameterizations, near quantitative agreement is found for SU(3) in D=2+l, but in D=3+l the 
extensions fail to remedy the qualitative disagreement. 



1. Introduction 

The traditional Isgur-Paton model Q was in- 
spired from the strong-coupling limit of lattice 
QCD. It describes a pure glue state as a closed 
loop of chromoelectric flux of vanishing diame- 
ter. The model is nonrelativistic and invokes an 
adiabatic limit such that motion is split between 
"fast" transverse oscillcations of the string and 
"slow" radial excitations. For SU{N > 3) the 
string carries an orientation vector. The predic- 
tions of the model depend greatly on the number 
of spatial dimensions, and so we will discuss the 
D=2+l case primarily. Adiabatically separating 
po allows quantization of the vibrational modes 
into phonons characterized by quantum occupa- 
tion number n^, where ± indicates the "helicity" 
of the mode. We also should expect a contribu- 
tion from the rest energy of the string, propor- 
tional to the string tension a. For a loop of radius 
PO) we have the Hamiltonian 

Tjr \ n I '^i,m''^i,m TTCT 2 2 /-, n 

H{pa) = 2Trcrpo + — + — 2_^m q,^^,{l) 

2TTapo 2po ^ 

where qi^m represents the normal coordinates for 
a,„ and (3m and TTi^m is the conjugate momen- 
tum. The zero-mode energy diverges but may 
be regularized by imposing a short-scale cutoff a 
on the string. Such a discretization leads to a 
zero-mode contribution as ^ — -rrr- + ■ • • ■ The 
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term linear in po can be identified with a " renor- 
malization" of the string tension a. The term 
linear in 1/po has been identified as a Liischer 
correction term, but inclusion of a rigidity term 
before quantization provides a l/pg dependence, 
and thus this term can be more easily identified 
with a renormalization of the effective elastic- 
ity in the extended model. The phonon modes 
can be shown to contribute to the energy as 
"^("■m + "-m)//°o, giving the mode energy factor 
^'^ = Hm^i'^tn + "-m)- For D=2+l, aU angu- 
lar momentum arises from the phonons, thus we 
have J = A = rn{n-^ - ). 

Since we expect the picture of the flux tube to 
break down near the origin, as well as to soften 
the singularity in the potential at po = 0, the po- 
tential in our effective Hamiltonian will be mul- 
tiplied by a heuristic suppression factor, F{p). 
Originally in 0, F{p) = (1 - e-f^P), where / 
is a free parameter near unity, but we find that a 
rational function, F{p) = (1 — //(p+ /)), gives a 
better fit to the lattice data. 

2. Inclusion of Elasticity 

Both intuition and calculation conclude that a 
flux tube should have nonvanishing thickness j^] 
in the transverse dimension. One manifestation 
of finite thickness we may expect to find even in 
our nonrelativistic model is an effective elastic- 
ity of the fiux tube. The free energy per unit 
length of a rigid bent rod is ^y^I, where y is 
the Young's modulus and / is the moment of in- 
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tertia perpendicular to the axis of bending. We 
expect y and / to be constant along the flux tube, 
thus integrating along the length of the flux tube 
f \yi-^dl cx ^ . In string language this is the 
curvature term. The constant of proportionality 
is our effective elasticity, denoted 7, and will be 
a free parameter in our action. We can see now 
that the effect of discretization on a scale a is 
to renormalize the elasticity 7 ^ 7 ~ jf • There 
has been some discussion as to the sign such 
an elasticity should take. We leave the elasticity 
as a free parameter letting the best fit parameter 
indentify the correct sign. 

3. Mixing 

Under charge conjugation, C, the loop orien- 
tation flips. To split the C = + and C = — 
sectors, we must introduce mixing between "left" 
and "right" orientated loops. We explore two 
scenarios: direct mixing and mixing through a 
higher topological state. Writing the indepen- 
dent orientations of the glue loop as t/^l and tpji^ 
the direct mixing matrix Hamiltonian is given in 
equation (||). For \I> = ( ^/j^, Vfl ) i 



H-L.R -a 



(2) 



where a is the mixing parameter and Hl^r is 
the Hamiltonian for a single glue loop. Since 
we would expect mixing to occur more easily for 
smaller loops, a may be multiplied by the comple- 
ment of the suppression function of the potential, 
so that the mixing goes as a(l — F{p)). Eigen- 
states for C are built as i/'l + V'fl and 1/;^ — ipR- 
We may instead include a state, denoted ipy, 
consisting of a glue loop with an extra diam- 
eter as an intermediate state through which 
one orientation may fluctuate to the other. 
Our matrix Hamiltonian then reads, for ^! — 



Hl,r 
Hl,r 
a rjca 



a 
■qca 
Hy 



(3) 



where the factor jyp = ±1 arises from coupling the 
ij^Y state to the appropriate admixture oi^L and 



ipR to match the C eigenvalue. Since the V'y state 
is not rotationally invariant, we must write the 
linear superposition as V'y = / rf6'e*'^^|'!/'Fe=o)- A. 
quick calculation shows that the il^y state has C 
eigenvalue Cy = (^l)"'- The combinations Vl ± 
ijjii have opposite C eigenvalues, thus by including 
the factor rjc = (—1)"' the ijiy state decouples 
from the admixture with the wrong C eigenvalue. 
As above, a may multiply a function of p, but the 
energetic cost of creating the diameter for large 
po provides an inherent suppression of mixing for 
large po- 

4. Building the Hamiltonian 

We build our Hamiltonian from (|l]) invoking an 
adiabatic separation of the variable p, as in the 
case of mesons [Q. When we take TTp^^ ~J^^ 
we identify our kinetic terms as 
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The potential F(po), given by the rest energy of 
the string and the contributions from the rigidity 
term and phonon modes, is 



V 
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The term 2v in is a contribution from the ver- 
tices each of static energy v. Such a contribution 
may follow from considerations in We also 
must include the rotational contribution for ipy, 



H, 



rotational 



2Iy 



J' 



{'iTT+±)api + 2vpl' 



(8) 



but wc do not include any phonon contributions. 

5. Results 

The spectrum is calculated by a suitable nu- 
merical diagonalization of the Hamiltonian and 
compared to recent lattice results (see, for exam- 
ple, 0). Wc minimize with respect to the 
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Figure 1 . Spectrum for direct mixing mechanism. 
/ = 0.17, 7 = 0.90, and a = 4.5. xVdof = 
44.1/(12-3) = 4.9. 
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Figure 2. Spectrum for Y state mixing mecha- 
nism. / = 0.28, 7 = 2.4, a = 3.8, and v = 9.3. 
xVdof= 64.7/(12-4) = 8.1. 



J = 0, 2 sectors by a simplex search in param- 
eter space and then calculate the full spectrum. 
Some ambiguity arises when identifying certain 
J^'-'' states for comparison. Lattice Monte Carlo 
simulation gives a below the l"*"^, whereas 
our model predicts the lightest ^ (at M = 8) 
well above the 1++. However, our model does 
predict a 4 below the 1++, suggesting that the 
Monte Carlo signal may be a J = 4 in disguise, 
a possibility for an Euclidean square lattice. Fig- 
ure |l| summarises our results for the direct mixing 
mechanism, with a multiplied by f / if + f) ■ Ex- 
cept for the 1^^, the model and the lattice are in 
remarkable agreement, although is somewhat 
larger than a believable fit should have. The de- 
generacy of the lattice data for the J = 1 sector 
is not available to the model, which maintains a 
significant splitting for all J such that the C — + 
states are lighter than the corresponding C — —. 

For mixing through Tpy, our results are shown 
in Figure |^, with no multiplicative suppression of 
a. Generally we find a poorer fit despite the ad- 
ditional parameter. Here we notice an interesting 
feature of this mechanism, namely that for even 
J the C = + states are lighter, but for odd J 
the C = — states are lighter. For the J = 1 sec- 
tor the lattice data must be averaged across par- 
ity to make our comparison, and the 1^^ might 
move relative to the 1=*=+ as better data become 
available. Which direction the 1^~ moves (if at 
all) would provide a reasonable criteria for distin- 
guishing between mixing mechanisms. 



6. Conclusions 

The Isgur-Paton model, extended by inclusion 
of an effective elasticity and by various mixing 
mechanisms, agrees reasonably with the lattice 
SU(3) spectrum in D=2-f 1, suggesting that a flux 
tube picture for glueballs might be viable there. 
In D=3-|-l, a qualitative discrepancy remains, as 
the model predicts a low-lying 1 ^ as an orbital 
excitation of the ground state 0++ which is at 
odds with the relatively heavy 1 ^ on the lat- 
tice. Nonetheless, the quality of the comparison 
in D=2-|-l suggests futher work be done mapping 
the pure gauge spectrum for this dimensionality. 
An interesting point for future comparison would 
be the mass of the J = 3 sector, which the model 
predicts to be signficantly lighter than the J = 1 
sector. 

We thank J. Paton for many useful discussions. 
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